Solar System tests do rule out 1/R gravity 
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Shortly after the addition of a 1/R term to the Einstein- Hilbert action was proposed as a solu- 
tion to the cosmic-acceleration puzzle, Chiba showed that such a theory violates Solar System tests 
of gravity. A flurry of recent papers have called Chiba's result into question. They argue that the 
spherically-symmetric vacuum spacetime in this theory is the Schwarzschild-de Sitter solution, mak- 
ing this theory consistent with Solar System tests. We point out that although the Schwarzschild-de 
Sitter solution exists in this theory, it is not the unique spherically-symmetric vacuum solution, and 
it is not the solution that describes the spacetime in the Solar System. The solution that correctly 
matches onto the stellar-interior solution differs from Schwarzschild-de Sitter in a way consistent 
with Chiba's claims. Thus, 1/7? gravity is ruled out by Solar System tests. 
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The discovery of an accelerated cosmic expansion |l|, 
2] has led to a flurry of theoretical activity. One class 
of solutions to the cosmic-acceleration puzzle consists of 
modifications to the general-relativistic theory of gravity. 
One particular proposal is the addition of a 1/R term to 
the Einstein-Hilbert action 0, Q . Such a term gives rise 
to a vacuum solution with constant curvature, the de 
Sitter spacetime, rather than the Minkowski vacuum of 
the usual Einstein-Hilbert action. 

Shortly after this proposal, Chiba []| argued that this 
theory is inconsistent with Solar System tests of gravity. 
In particular, he showed that the theory is equivalent to a 
scalar-tensor theory that is known to make Solar System 
predictions that conflict with measurements. 

Since then, however, there have been a number of pa- 
pers arguing or implying that Chiba's analysis is flawed 
H El 03 ■ The crux of the counter-argument is that 
1/i? theories admit as a static spherically-symmetric so- 
lution the usual vacuum Schwarzschild-de Sitter space- 
time. Apart from a cosmological constant that is too 
small by many orders of magnitude to affect anything ob- 
servable in the Solar System, these solutions are just the 
usual Schwarzschild solution. Consequently, they argue, 
there is no effective difference between the Solar System 
spacetime in these models and that in ordinary general 
relativity. 

Here we point out that these arguments are incor- 
rect, and that Chiba was right. The crucial point is 
that although the Schwarzschild-de Sitter spacetime is 
indeed a spherically-symmetric vacuum solution to the 
1/R equations of motion, it is not the unique spherically- 
symmetric vacuum solution in this theory. The correct 
solution is determined by matching onto the solution in 
the interior of the star. When this is done correctly, 
it is found that the Schwarzschild-de Sitter spacetime 
does not describe the spacetime around the Sun, and 
that Chiba's result stands. This misunderstanding has 
now propagated through a number of papers. There are 
moreover a number of other papers that cite these incor- 
rect papers in a way that suggests that they may be onto 



something. We thought it worthwhile to correct the error 
before it propagates any further. 

Before describing the correct spherically-symmetric 
spacetime for 1/R gravity, we consider a very simple 
and analogous problem that illustrates what is going on. 
Suppose we wanted to know the electric field around a 
spherically-symmetric charge distribution p{r) confined 
to radii r < R. For radii r > R, the Poisson equation 
V 2 = 47T/9 relating the electric potential (j> to the charge- 
density distribution p reduces to V 2 </> = 0. A spherically 
symmetric solution to this equation, one might argue, is 
4> = 0, implying no electric field. This is clearly incorrect. 

What went wrong? Although <f> = is indeed a spheri- 
cally symmetric solution to V 2 = 0, it is not the unique 
solution. Another solution is <j> = c / r i f° r r > R. The 
constant c in this equation is furthermore fixed in this 
case to be c = Q, where Q = J p d 3 x is the total charge, 
by integrating the right- and left-hand sides of the Pois- 
son equation V 2 </> = Airp over the entire volume. 

In brief, something similar happens in 1/R gravity. 
The differential equations for the metric components 
gu(r) and g rr {r) are supplemented by a differential equa- 
tion for the curvature R, as we will see below. The three 
differential equations have the Schwarzschild-de Sitter 
spacetime as a solution, but these vacuum solutions do 
not match onto the solutions in the presence of a source 
(i.e., the Sun). There is an additional vacuum solution 
that correctly matches onto the solution in the presence 
of the source. 

Now the details: the gravitational action of 1/R grav- 
ity, 

(i) 

may be varied with respect to the metric g M „ to obtain 
the field equation Q 

^GT, V = (l + ^) *W - \ (l - |i) R 9»» 

+/J 4 (.g^V«V Q - V^V,) R- 2 . (2) 
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We begin by using the trace of the field equation to de- 
termine the Ricci scalar R. Contracting Eq. J5J with the 
inverse metric yields 
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where T = g^T^. 

The constant-curvature vacuum solution is obtained 
by setting T = and V ^R = 0. It is R 2 = 3/i 4 , 
corresponding to the de Sitter spacetime with Hubble 
parameter H 2 = /i 2 /(4-\/3), equivalent to the general- 
relativistic vacuum solution with a cosmological constant 
A = 3H 2 = \/3p 2 /4. The metric for this spacetime can 
be written as a static spherically-symmetric spacetime: 



ds 2 = - (1 - H 2 r 2 ) dt 2 + (l - H 



' dr 2 +r 2 dn 2 . (4) 



To match the observed acceleration of the universe, the 
effective cosmological constant must be set to A ~ p 2 ~ 
H 2 ~ lO" 56 cm" 2 . 

We now consider the spacetime in the Solar System 
in this theory. First of all, the distances (~ 10 13 cm) 
in the Solar System are tiny compared with the distance 
p^ 1 ~ 10 28 cm, so fir <C 1 everywhere in the Solar Sys- 
tem. Moreover, the densities and velocities in the Solar 
System are sufficiently small that we can treat the space- 
time as a small perturbation to the de Sitter spacetime. 
The spacetime should also be spherically symmetric and 
static. The most general static spherically-symmetric 
perturbation to the vacuum de Sitter spacetime given 
by Eq. (0 can be written 

ds 2 = - [1 + a{r) - H 2 r 2 ] dt 2 

+ [1 + b(r) - H 2 r 2 Y X dr 2 +r 2 dtt 2 , (5) 

where the metric-perturbation variables a(r),b(r) <C 1. 
In the following, we work to linear order in a and b, and 
also recall that fir <C 1. However, a, b are not necessarily 
small compared with fir. 

We now return to the trace of the field equation, given 
by Eq. (0J), and solve it for the Ricci scalar R(r) in the 
presence of the Sun. We write the trace equation in terms 
of a new function, 
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c(r) = + 
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R 2 (r) : 



(6) 



and demand that c(r) — ► as r — > oo so that R ap- 
proaches its background value of \/?>p 2 far from the 
source of the perturbation. Therefore, c(r) parameter- 
izes the departure of R from the vacuum solution, and 
we anticipate that c(r) will be the same order in the per- 
turbation amplitude as the metric perturbations a(r) and 
b(r). In terms of c(r), Eq. @ becomes an exact equation, 
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In the Newtonian limit appropriate for the Solar Sys- 
tem, the pressure p is negligible compared to the energy 
density p, and so T = —p. Neglecting terms that are 
higher order in a(r), b(r), and p 2 r 2 , we are able to rewrite 
Eq. 10 as 



V 2 c + V3p 2 c = — p, 



(8) 



where V 2 is the fiat-space Laplacian operator. Note that 
in writing this equation, which is linear in c(r), we have 
also neglected higher-order terms in c(r). Below, we will 
check that the solutions we obtain have c(r) <C 1 ev- 
erywhere, consistent with our assumptions. The Green's 
function for Eq. (JSJ is — cos(3 1 / 4 /Ltr)/(47rr). Convolving 
this with the density gives us the solution to Eq. |JSJ. 
However, we are restricting our attention to the region 
where pr <C 1 , so the Green's function reduces to that for 
the Laplacian operator. Therefore the equation we need 
to solve is V 2 c = — (8nGp)/3. Integrating the right- 
hand side over a spherical volume of radius r gives us 
— 8nGm(r)/3, where m(r) is the mass enclosed by a ra- 
dius r. Using Gauss's law to integrate the left-hand side 
gives us 47rr 2 c'(r), where the prime denotes differentia- 
tion with respect to r. Thus, the equation for c(r) be- 
comes 

dc _ 2Gm(r) 
dr 3r 2 

Integrating Eq. 10 and using the boundary condition 
that c — > as r — > oo gives us the solution c(r) = 
(2/3)(GM/r)[l + 0{pr)} for r > R Q . Note also that 
integration of the equation for c'(r) to radii r < Rq in- 
side the star implies that the scalar curvature R remains 
of order p 2 , even inside the star. We thus see that c <C 1, 
so we were justified in using the linearized equation for 
c(r). 

This solution for c(r) implies that 
GM^ 
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R = V3p 2 
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We have thus shown that R is not constant outside the 
star and have already arrived at a result at odds with the 
constant-curvature Schwarzschild-de Sitter solution. No- 
tice that had we {incorrectly) used p = in Eq. |8|l. then 
the equations would have admitted the solution c(r) = 0; 
i.e. the constant-curvature solution. However, this would 
be incorrect, because even though p = at r > i?©, the 
solution to the differential equation at r > Rq depends 
on the mass distribution p(r) at r < i? Q . In other words, 
although the Schwarzschild-de Sitter solution is a static 
spherically-symmetric solution to the vacuum Einstein 
equations, it is not the solution that correctly matches 
onto the solution inside the star. Note further that the 
solution for R both inside and outside the star is (to lin- 
ear order in c), 



R = V3p 2 



1 - -dr) 
2 y ' 



(11) 
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Clearly, 1/ R gravity produces a spacetime inside the star 
that is very different from general relativity. This result 
shows that in this theory one should not assume that 
R = 8irGp; this has lead to some confusion 0, 0, . 

To proceed to the solutions for a(r) and b(r), we re- 
arrange the field equation for 1/R gravity [Eq. J2J] to 
obtain equations, 



R 2 
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for the Ricci tensor in terms of the Ricci scalar. When 
the expression for R obtained from the trace equation is 
inserted into the right-hand side, we obtain equations for 
the nonzero components of the Ricci tensor, 

|v 2 c, (13) 
(14) 

- + c"(r)\ (15) 



where we have neglected terms of order /i 2 c, Gpc and c 2 
in all three expressions. 

For the perturbed metric given by Eq. 10), the tt 
component of the Ricci tensor is (to linear order in 
small quantities) R\ = 3H 2 - (l/2)V 2 a(r). Applying 
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tion for a(r), 



inGp)/3 to Eq. (|13|) leaves us with an equa- 
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(16) 



plus terms that are higher order in GM/r and pr. The 
solution to this equation parallels that for c(r); it is 



(17) 



dr r* 

both inside and outside the star. Outside the star, this 
expression may be integrated, subject to the boundary 
condition a(r) — > as r — > oo, to obtain the metric per- 
turbation, 

2GM , , 

air) = , r > R Q , (18) 

r 

exterior to the star. Note that this recovers the Newto- 
nian limit for the motion of nonrelativistic bodies in the 
Solar System, as it should. 

The rr component of the Ricci tensor is (to linear order 
in small quantities) R r r = 3H 2 - (b'/r) - (a"/2). Given 
our solution for a'(r) and c'(r) = — (2/3)Gm(r)/r 2 , 
Eq. i|14|) becomes a simple differential equation for b(r), 

db Gm(r) Gm'(r) 
dr r 2 r 



d_ 

dr 



-Gm{r) 



(19) 



Integrating this equation subject to the boundary condi- 
tion b{r) — > as r — > oo gives an expression for b(r) that 
is applicable both inside and outside the star: 



6(r) 



Gm(r) 
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This expression for b{r) and Eq. (|17|) for a'{r) also satisfy 
Eq. H15|) for the angular components of the Ricci tensor. 
The Ricci scalar [Eq. (|1L)|) ] is recovered from the Ricci 
tensor components if terms higher order in 0(pr 2 GM/r) 
are included in our expressions for a(r) and b(r). 
The linearized metric outside the star thus becomes 



r 

GM 



(21) 



H 2 r 2 ) dr 2 + r 2 dQ 2 . (22) 



Noting that in the Solar System, Hr ^ 1 and that the 
PPN parameter 7 is defined by the metric, 



ds 2 = - 1 - 



2GM^ 



dt 2 



l + 2 -^-)dr 2 + r 2 dn 2 , 



(23) 

we find that 7 = 1/2 for 1/R gravity, in agreement 
with Chiba's claims H,0]> and prior calculations; e.g., 
Refs. 0, 0. We note that recent measurements give 
7=1 + (2.1 ± 2.3) x 1CT 5 003. 

Other authors have noted that Birkhoff's theorem — 
that the unique static spherically-symmetric vacuum 
spacetime in general relativity is the Schwarzschild 
spacetime — is lost in 1/R gravity, and that there may 
be several spherically-symmetric vacuum spacetimes. Al- 
though this is true, what we have shown here is that the 
Solar System spacetime is determined uniquely by match- 
ing the exterior vacuum solution to the interior solution. 
When this is done correctly, it is found that the theory 
predicts a PPN parameter 7 = 1/2 in gross violation of 
the measurements, which require 7 to be extremely close 
to unity. 

A few final comments: It is important to note that 
the structure of 1/R gravity (for example, the way mat- 
ter sources the metric) is completely different than the 
structure of general relativity, even in the limit p — > 0. 
In particular, the theory does not reduce to general rela- 
tivity in the p — > limit, and this can lead to confusion. 
This is due to the fact that the introduction of additional 
terms in the Einstein-Hilbert action brings to life a scalar 
degree of freedom that lies dormant in general relativity. 
We also note that Chiba's mapping of f(R) theories to 
scalar-tensor theories is perfectly valid; it amounts to no 
more than a variable change, from R to <f> = 1 + p 4 /R 2 . 
The trace equation, Eq. ©, is then equivalent to the 
scalar-field equation of motion in the scalar-tensor the- 
ory. Also, the fact that general relativity is not recov- 
ered in the p — » limit becomes particularly apparent in 
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the scalar-tensor theory, as we will discuss elsewhere. Al- 
though we have restricted our analysis, for clarity, to 1/ R 
theory, similar results can also be derived for other f(R) 
theories. For example, the correct matching of the exte- 
rior and interior solutions can be used to distinguish be- 
tween the spherically-symmetric vacuum spacetimes for 
R 1+s gravity discussed in Ref. 0. We plan to present 
more details in a forthcoming publication [2fj | . 
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